Global phase diagram and six-state clock universality behavior in the triangular 
antiferromagnetic Ising model with anisotropic next-nearest-neighbor coupling: 

Level-spectroscopy approach 
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We investigate the triangular-lattice antiferromagnetic Ising model with a spatially anisotropic 
next-nearest-neighbor ferromagnetic coupling, which was first discussed by Kitatani and Oguchi. By 
employing the effective geometric factor, we analyze the scaling dimensions of the operators around 
the Berezinskii-Kosterlitz-Thouless (BKT) transition lines, and determine the global phase diagram. 
Our numerical data exhibit that two types of BKT-transition lines separate the intermediate critical 
region from the ordered and disordered phases, and they do not merge into a single curve in the 
antiferromagnetic region. We also estimate the central charge and perform some consistency checks 
among scaling dimensions in order to provide the evidence of the six-state clock universality. Further, 
we provide an analysis of the shapes of boundaries based on the crossover argument. 
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I. INTRODUCTION 

The antiferromagnetic Ising model on the triangular 
lattice may be the simplest example to possess the frus- 
tration effects. As its manifestation, the system does not 
complete the long-range order even in the ground state, 
but has the critical ground-state ensemble HI 0, IS El El • 
This emerging criticality due to the frustration is the 
Gaussian type, where the Coulomb-gas picture @ well 
describes properties of elementary excitations and also 
singularities of physical quantities Meanwhile, 
continuum field theories for systems nearby the criticality 
have been well developed on the basis of the conformal 
symmetry; it offers, beyond the Coulomb-gas descrip- 
tion, the unified and pervasive approaches to investigate 
the one-dimensional (ID) quantum and two-dimensional 
(2D) classical systems. Further, since they also offer us 
powerful strategies to numerically investigate the systems 
around the Gaussian fixed point, quantitative reliability 
becomes higher than ever before, which contributes to 
clarify lon^-sta nding problems in related research fields 

In this paper, we treat the triangular antiferromagnetic 
Ising model (TAFIM) with anisotropic next-nearest- 
neighbor (NNN) ferromagnetic coupling which wc 
call the Kitatani-Oguchi model hereafter. The reduced 
Hamiltonian including temperature TL = 0H is given by 



H{K X ,K 2 ) 



K, 



m' 



(i) 



The binary variable o~j = 0, 1 is on the jth site of the tri- 
angular lattice A which consists of interpenetrating three 
sublattices A; (/ = 0,1,2). The first sum runs over all 
nearest- neighbor (NN) pairs (J, k) , whereas the second 
over NNN pairs [j, k]' in two of the three directions (see 
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FIG. 1: The schematic representation of Hamiltonian Q. 
The jth site is specified by two integers (ji, J2) as labeled 
in the figure. The dotted (dashed) lines show the NN AF 
(anisotropic NNN F) coupling. An example of the spin con- 
figurations with the y/3 x \/3 structure is exhibited, where the 
spins on two (one) of three sublattices Ao,i (A2) are up © 
(down Cg>). 



Fig. Although the model was introduced due to a 
technical reason in their numerical calculation, unlike in 
the isotropic NNN F coupling case, it provides an oppor- 
tunity to investigate an important effect, i.e., the spatial 
anisotropy in critical phenomena. 

Since its introduction, there is a belief that the model 
is not strongly anisotropic and belongs to the six- 
state clock (6SC) universality class as the isotropic model 
does (for more severe anisotropies, see 0,^3). Thus, 
two types of Bcrezinskii-Kostcrlitz-Thouless (BKT) tran- 
sition lines are expected to separate the critical region 
from the ordered and disordered phases [l8L Il9l | . Actu- 
ally, Kitatani and Oguchi (KO) performed the transfer- 
matrix (TM) calculations at the ratio K2/K1 = 0.5 in 
order to evaluate the function within the Roomany- 
Wyld approximation |20j and the spin-spin correlation 
function, and then concluded the existence of the critical 
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finite temperature region Subsequently, Miyashita, 
Kitatani, and Kanada performed large scale Monte Carlo 
(MC) simulation calculations and found the consistent re- 
sults |2l| . Using the cluster TM method (an approxima- 
tion), Pajersky and Surda stated that the intermediate 
phase possessed an incommensurate nature |22|. Con- 
trary to these, de Queiroz and Domany performed the 
TM calculations, and reported limited evidence for an ex- 
istence of the BKT phase on the basis of the phenomcno- 
logical renormalization-group (PRG) analysis j2^|. They 
argued that this limitation was due to the lack of sixfold 
symmetry, and they also provided the qualitative phase 
diagram with a multicritical point. Quite recently, Qian 
and Blote succeeded in obtaining the phase diagram of 
TAFIM with isotropic NNN F coupling using both TM 
and MC calculations, whereas, for present model, they 
stated that an app lication of the conformal mapping be- 
came difficult |24|. Considering all these together, there 
still exist some unclear points due to the nature of the 
BKT transitions and/or the effects of the anisotropy. 

In this paper, we shall provide the global phase dia- 
gram of the Kitatani-Oguchi model Q on the basis of 
the TM calculation data. Note especially, our result ex- 
hibits that the phase diagram possesses the same struc- 
ture as that of the isotropic model 13- Therefore, 
this provides evidence to support KO's motivation and 
confirms the present model being in the 6SC universality 
class, independently of the spatial anisotropy. 

The organization of this paper is as follows: In Sec. 
ITT1 we shall discuss the effective field theory to describe 
the low-energy and long-distance behaviors of the model, 
where the effects of the NN and the anisotropic NNN 
couplings are explained. Simultaneously, since we shall 
employ the level-spectros copy method to treat the two 
types of BKT transitions [25j , we address some relevant 
issues necessary for its application. In Sec. IIIII we shall 
give our numerical estimates of the phase boundary lines. 
By introducing the parameter-dependent geometric fac- 
tor to obtain the isotropic description of the model in 
the 2D Euclidean space and by evaluating the factor on 
the basis of the conformal field theory (CFT), we check 
the criticality of the system and some universal relations 
among the scaling dimensions. The last section, Sec. II VI 
is devoted to the discussion and the summary of our re- 
search. We provide the analysis of the shapes of bound- 
aries based on the crossover argument. We also compare 
our data with previous research results, and then give 
some comments. 



II. THEORETICAL DESCRIPTION AND 
NUMERICAL CALCULATION METHOD 

A. Effective field theory in the sine-Gordon 
language 

We shall start by referring to some relating results in 
the literature. In the exactly solved case K2 = 0, the 



system shows the 2D Gaussian criticality at K\ = 00. In 
the scaling limit of the lattice spacing a — > 0, but keep- 
ing x = (£1,3:2) = (aji, aj2/C) finite [£ = 2/V3 is the 
geometric factor for the triangular lattice], its effective 
description is given by the Lagrangian density 
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where K (= i) is the Gaussian coupling and the phase 

field satisfies the periodicity condition \ 2(f) + 2tt = \/2<p 
El ■ Reflecting the discreteness of height variables hj 
in the triangular Ising solid-on-solid model which is ob- 
tained from the zero-temperature TAFIM , there exists 
the nonlinear potential cos6\/2(f>. However, it is highly 
irrelevant at the point 0] , so here we dropped it (see 
below) . Writing Sj = e l7r<Tj and the distance between jth 
and fcth sites rjk = |x — x'|, then the asymptotic behavior 
of the spin-spin correlation function is given as {sjSk) — 

A cos ip^/r 1 ^ + B/r®[ 2 , where ip jk =ipj — ip k and a sub- 
lattice dependent phase ifj := J2i=o 1 2 SieAj (^V^^ij 
(A, B constants.) 0,0. Thus, the scaling dimensions of 
the uniform and staggered magnetizations s and S con- 
structed through a coarse graining of Sj and e ±%Vl Sj with 
respect to each elementary triangle are given by x s = | 
and xs = j, respectively. 

Next, we move on to the anisotropic case A'2 7^ 0. 
As asserted by KO, there is a plausible reason to expect 
that anisotropic NNN F coupling plays the same role as 
the isotropic one 0|, i.e., it only reduces the Gaussian 
coupling down to a certain value 0. However, there is 
no reason to expect that £ = 2/ V3 is appropriate to 
obtain a rotationally invariant description of the model 
in the 2D Euclidean space. This is because £ = 2/V3 
can express the geometric structure of the lattice A, but 
the anisotropic coupling requires a further rescaling fac- 
tor, which must be incorporated multiplicatively into the 
geometric factor. Therefore, we suppose that the effec- 
tive value of the geometric factor £ is not known a pri- 
ori and depending on the coupling values. In some sys- 
tems including those in somewhat different situation, the 
anisotropy effect has been discussed 0, ^| , and fur- 
ther it is also important in the ID quantum systems de- 
scribed by the Tomonaga-Luttinger-liquid Hamiltonian 

M 



H Thh (v,K) 



dx- 



2tt 



k (d x ef + k- 1 (d x $) 



(3) 



where [$>(x),d X '0{x')/Tv] = i8(x — x'), and v is the veloc- 
ity of an elementary excitation. While, in the imaginary- 
time formalism, the (1 + 1) space-time is treated as 2D 
Euclidean, x = (x,vt), so as to extract the Lagrangian 
density @|, v as well as K depends on the lattice- model 
parameters. Returning to our case, as we will see be- 
low, this dependence indeed requires an extra procedure 
to estimate £ in numerical calculations, but this changes 
nothing of the theoretical description, provided £ is cho- 
sen properly. Therefore, we advance our argument. 
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For large enough K2, the Gaussian criticality is lifted 
and the sixfold-degenerate ordered state with the y/3xV3 
sublattice structure is stabilized (see Fig. P). On the 
other hand, for small K2, the thermal scaling field u = 
e -Ai introduces the relevant energy-density perturbation 
(its bosonized expression is e = y/2 cos y/29 in terms of 8 
being dual to 4>) and brings about the disordered state • 
In summary, these are described by the dual sine-Gordon 
model whose Lagrangian density is given by 

C = C + -^—^ (y+ cos 6v^> + y g cos V26») , (4) 

where y$ < (see below) and ye <x u are the coupling 
constants |27| ■ a is a short-distance cutoff. Since the 
nonlinear terms in Eq. 10} are both irrelevant in the re- 
gion | > K > the present model (|TJ may possess 
two BKT-transition lines We define v = e K2 and 

denote the upper- and lower-temperature boundaries as 
v\j(u) and vi_,(u) [note v\j(u) < i>l(u)] for convenience, 
then their precise determinations arc our main goal. 

Quite recently, Matsuo and Nomura investigated the 
classical 2D 6SC model 0. Especially, mapping it to 
the ID quantum 6SC model with the explicit duality re- 
lation, they succeeded in determining two BKT transi- 
tion points, where the crossings of levels observed in the 
systems with periodic and twisted boundary conditions 
were used. Although their criteria may also be efficient 
to our case, an implementation of the twisted boundary 
condition is unclear at present. On the other hand, we 
have also treated the BKT transitions observed in the 
AF three-state Potts model with the NNN F coupling by 
employing rather naive criteria Therefore, we shall 
employ the same approach in the following (see also H3)- 

B. Finite-size estimates of BKT-transition points 
in upper and lower temperatures 



the BKT-transition point t = y 1 /y - 1. On the BKT 
line, yo = yi — On the other hand, another 

important operator is the uniform magnetization, s — 
y/2 cos 3y/2(f> in its bosonized form, whose dimension is 
x s (l) ~ to (2 — yo) in the same region. Consequently, 
the level-crossing condition, 

xo{l) = yZ s (Z), (6) 

provides a finite-size estimate of v\j(u), v\j(u,L). 

Next, we consider a region near vt,(u) where cos\f20 
is irrelevant. Then, the effective Lagrangian density is 

In this case, M and y/2 cos 6y/2(f> hybridize and result in 
the Al-like and cos-like operators, O2 and O3, respec- 
tively. By redefining the coupling constants {yoiVi) = 
(18 K — 2, — y^), we obtain the scaling dimensions X2(l) — 
2 — Vo (l + §*) an d x 3 (l) ~ 2 + ?/o (2 + |t) near the tran- 
sition point. Since the uniform magnetization has the 
dimension x s (l) ~ \ [2 — yo (1 + 2i)], the level-crossing 
condition for vi,(u,L) may be given by 

x 2 (l) = Ax s (l). (8) 

Here, let the following be observed: Since we have sup- 
posed y$ < 0, the sixfold-degenerate ordered states cor- 
respond to the six locking points of the phase variable, 
i.e., (\f24>) ~ nq/3,(q = 0, 1, ■ ■ ■ ,5). According to the 
bosonized expressions, the averages of the uniform and 
staggered magnetizations [S = exp(±iv2<^)] take values 
(s) cx e %lxq and (S) oc e ±47ri?/3 in the ordered states. These 
are consistent with the relationship between the phase 
values and the spin configurations given in Refs. [71. l28j. 



First, let the system around v\j(u) be considered where 
cos 6V20 is irrelevant. Then, it is well described by 

With respect to the determination of the BKT-transition 
points, one of the authors (K.N.) pointed out the im- 
portance of marginal operators |25| : Especially, in this 
case, M = (1/K) (V0) 2 and e = y/2 cosy/29 hybridize 
along the RG flow and result in two orthogonalized op- 
erators, i.e., the ".M-like" and the "cos-like" operators. 
We denote the former and latter as Oq and 0\, and de- 
fine the system on A with M (— > 00) rows of L (a mul- 
tiple of 3) sites wrapped on a cylinder, say A(L x 00). 
Then, according to the conformal perturbation theory, 
their renormalized scaling dimensions near the multicrit- 
ical point (y ,yx) = (1/2K - 2, ye) = (0,0) are given 
as x (l) ~ 2 - y (l + ft) and Xl {l) ~ 2 + y Q (2 + |f) , 
respectively, where I = InL and a small deviation from 



C. Discrete symmetry properties of excitations 

Now, let us consider the system on A(L xoo), and de- 
fine the transfer matrix, T(L), connecting the NNN rows 
in the vertical direction of Fig. 2] We denote its eigenval- 
ues as X P (L) or their logarithms as E p (L) = —5 In |A p (L)| 
(p specifies a level) . Then, the conformal invariance pro- 
vides direct expressions of the central charge c, the scal- 
ing dimension x p and the conformal spin s p in the critical 
systems as follows [53, 

E S (L) ~ Lf - g^c, (9) 

AEp(L) ~ ^x p , Ak p (L) ~ ^-s p . (10) 

Here, E S (L), AE p (L) [= E p (L) - E S (L)], and Ak p {L) 
correspond to the ground-state energy, an excitation gap, 
and a momentum, respectively. / is free energy per site. 
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Furthermore, the following formulas are available for the 
Gaussian system: 



1 



x 



Kn z 



nr 



V"" ' K 
Sp = nm+(N-N), 



+ (N + N), (11) 
(12) 



where (n, m) are the electric and the magnetic charges 
@ , and (iV, N) are non-negative integers [3ll ■ 

In the numerical calculations of the above-mentioned 
excitation levels, the symmetry properties such as the 
translation of one lattice spacing (a cyclic permutation 
among sublattices A 4 ) T, the space inversion V , and the 
spin reversal S are quite important. This is because these 
symmetry operations can be also interpreted in the field 
language as T : V2<f> i-> y/2(j> + 2_7r/3, V : \f2<f) h-> — \f2(f), 
and S : \/2c/> i-> V2<f> + tt Q, U |2g. Therefore, the cor- 
responding levels to marginal operators Oq and Oi (O2 
and O3) can be found in the subspace of the wavenum- 
bcr k = and the even parity for both V and S. On the 
other hand, the uniform magnetization s = \/2 cos 3 \2<f) 
is k = for T and even for V, but it is odd for S as ex- 
pected. We thus calculate the excitation levels AE P (L) 
by utilizing these symmetry operations and solve the 
level-crossing conditions © and (jSJ), numerically. 

Here, the following should be remarked on. When us- 
ing the KT criterion to determine the BKT-transition 
points [e.g., K = j for Vjj(u)], we should estimate the 
Gaussian coupling from an appropriate excitation gap 
through the former of Eq. (|10|) . If, like the present case, 
C is not known a priori, this requires its estimate in ad- 
vance of the gap data. On the other hand, since the 
level-crossing conditions 10 and (jSJ) arc homogeneous ex- 
pressions in terms of the scaling dimensions, the corre- 
sponding excitation gaps (or scaled ones) , instead of the 
scaling dimensions, can be used to estimate the BKT- 
transition points. This property is one of the advantages 
of the level-spectroscopy approach in the studies of ID 
quantum and anisotropic 2D classical systems. 
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FIG. 2: The K2 dependence of the scaled gaps at K\ — 1 (the 
correspondences between marks and system sizes are given in 
these panels), (a) [(c)] shows Xo and 16X s /9 {Xo 2 and 
AXs). Crossing points give finite-size estimates vu,l(m, L). 
(b) gives data Xm and 8X s /3 whose crossing points estimate 
the self-dual line v B d(u,L). 
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III. RESULTS 

A. Phase boundaries and the self-dual line 

Now we perform the exact-diagonalization calculations 
of T(L) for systems up to L = 24 by the use of the 
Lanczos algorithm. In Figs. |2Ja) and Hfc), we pl°t 
examples of the K2 dependencies of the scaled gaps 
X p (L) :— AEp(L) / (2ir / L) (or values multiplied by con- 
stants for convenience) at K\ = 1. Then, we can find the 
points at which the above condition © or © is satisfied 
[i.e., Utj,l(w, L)]. In addition to the logarithmic correc- 
tions, there is another type of correction stemming from 
the x = 4 irrelevant operators [32l | , we shall thus extrap- 
olate them to the thermodynamic limit according to the 
least-squares fitting of the polynomial in 1/L 2 . 



FIG. 3: The phase diagram. Open circles, squares and tri- 
angles with curves exhibit the BKT-transition lines vxj(u), 
vl(u), and the second-order transition line Vi(u), respectively. 
The thick vertical line gives the first-order transition bound- 
ary. The filled circle at (1,1 + \/2) denotes the decoupling 
point with the three independent Ising criticality. Diamonds 
with the dash-dotted line shows v s d(u). 



In Fig. |31 we give our phase diagram in the 2D model 
parameter space (u,v) = (e~ Kl , e K ' 2 ). The open circles 
and squares with the solid curves exhibit the lines v\j(u) 
and i>l(m), respectively, and they separate an intermedi- 
ate region from the disordered phase and from the or- 
dered phase "Ordered (I)" with sixfold degeneracy. The 
filled circle at (u, v) = (1,1 + V%) denotes the decoupling 
point with three independent Ising criticality. To com- 
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(b) 



TABLE I: Examples of the L dependences of the scaling 
dimensions x\ and xa and the averages x av and x' av (see the 
text) on the BKT-transition points i>u,l(u) {K\ — 1). We 
extrapolate the finite-size estimates to L — > oo using the least- 
squares fitting of the polynomial in 1/L 2 . 



L 



12 



15 



18 



21 



24 



zi(/) 2.48863 2.50955 2.51917 2.52263 2.52333 

■gavffl 1.87339 1.90402 1.92525 1.94038 1.95161 1.988 

x 4 (l) 0.83533 0.80457 0.78318 0.76714 0.75449 

Xa V (0 0.49797 0.49447 0.49202 0.49014 0.48862 0.484 



0.9 ■ 1 ■ 1 ■ 1 ■ 1 

05 (v-1)/V2 1 

FIG. 4: The inverse effective geometric factor £ _1 (the central 
charge c) versus (v — l)/V% is plotted in the upper panel (a) 
[the lower panel (b)]. Filled marks plot the data with u — 0, 
and open ones along v E d(u). £ _1 takes the close value to \/3/2 
at v = 1. 

plete the phase diagram, we also calculate the second- 
order phase transition point V\[u) with the Ising critical- 
ity in the F region u > 1, where the finite-size estimates 
by the PRG method are extrapolated to the thermody- 
namic limit according to v\(u,L) ~ vjju) + a/L 3 (see 
the triangles with the solid curve) |U |^ ■ The ordered 
phase "Ordered (II)" has the twofold degeneracy, and the 
thick vertical line shows the first-order phase transition 
boundary between ordered phases. 

For a later discussion, here we mention the self-dual 
line v s d(u) embedded in the critical region (3^. Al- 
though it is not the phase transition lines, it is expected 
to be good for numerical calculations [T3 |. Defining the 
transformation 60 «-> 6, then we can find the duality 
relation of the effective model (@J), i.e., {K.y^.yo) «-> 
(1/3QK, yg, y^,). Thus it becomes invariant at K = i and 
V<P — ye- Since the self duality provides the degeneracy of 
excitation levels, e.g., cos6^/2(j) an d cosv2#, their cross- 
ing provides the finite-size estimates v s d(u,L). However, 
these are higher energy excitation levels, and their stable 
estimations by the Lanczos method are rather difficult. 
Alternatively, we employ the condition x_m(1) = 8x s (l)/3 
[%m{ = 2) is the dimension of Fig. E^b) exemplifies 
the level crossing. Extrapolating them to the thermody- 
namic limit as v s d(u, L) ~ v s d(u) + a/L 2 , we determine 
v s d(u) (diamonds with the dash-dotted line in Fig. |2J, 
which is between two boundaries wtj,l(w), and is termi- 
nated at the critical decoupling point. 



B. Effective geometric factor, central charge, and 
consistency checks 

Although our calculations so far do not need the nu- 
merical estimation of it is necessary to check the crit- 
icality and the universal relations among scaling dimen- 
sions. For this, following the recent development in the 



study of ID quantum systems |36j, we focus our atten- 
tion to the level- 1 descendant in the conformal tower of 
the identity operator, i.e., in the CFT language. 

Since the corresponding excitation level is specified by 
(n, m, N, N) = (0, 0, 1, 0), we can estimate £ as 

AE f JL,k = 2ir/L) 
C 1 = lim — . (13) 

L^ca 2TT I L 

The finite-size estimates, i.e., the right-hand side ratio, 
are extrapolated to the thermodynamic limit by the least- 
squares fitting of the polynomial in 1/L 2 . Then, we ob- 
tain the results which arc plotted in Fig.^Ja) (the upper 
panel). Filled triangles plot the data along the u = 
line, and open ones along v s a(u). Here note that since 
Eq. I|13|) is only valid for the system with the Gaussian 
criticality, it does not hold for others including the crit- 
ical decoupling point. For the isotropic case v = 1, the 
estimation excellently agrees with the exact value, and 
it increases with v as expected. On the other hand, for 
three decoupled square lattices formed by the anisotropic 
NNN coupling, C _1 may equal to |, so that £ seemingly 
jumps at the critical decoupling point. This might have 
relevance with the jump of the central charge, but a more 
detailed analysis is left for future study. 

Now, according to Eq. 10, we can estimate the central 
charge c from the L dependence of the ground-state en- 
ergy E g (L); the results are plotted in Fig.^Jb) (the lower 
panel). Although £ -1 increases to nearly twice as large 
as the isotropic value, the central charge keeps c = 1 
within 1.5% deviations, which clearly demonstrates the 
Gaussian criticality of the model. 

Next, we shall check some relations. Since the am- 
plitudes of the logarithmic corrections are given by the 
operator-product-expansion coefficients, some universal 
relations among the scaling dimensions have been dis- 
covered: For instance |25ll37j. 

i[2x (0+Zi(0]-2 onwuH, (14) 
\[3x s {l) + x 4 {l)} ~ | onv L (u). (15) 

Here, x^il) is the scaling dimension of v^sinSv^t/), which 
is k = for T and odd for V and S. In Table we 
give, as an example, the dimensions at Ki = 1 esti- 
mated using the former of Eq. (|lUp . Although Xq(1) and 
x\{l) [x s {l) and xn(l)] largely deviate from the free-boson 
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value 2 (i) due to the logarithmic corrections, their main 
parts cancel each other in Eqs. (|14fl and 1)150. Therefore, 
the average x av (x' av ), the left-hand side of Eq. ifFf]! [Eq. 
(|15fl ], takes the value close to 2 (i). These checks can 
be passed only if the systems are on the BKT-transition 
lines vu,l(m), and the numerically utilized levels have the 
theoretically expected interpretations. Therefore, these 
arc helpful to demonstrate the reliability of our numerical 
results. 



C. Summary of results 

Consequently, we can confirm that the intermediate 
phase shows the Gaussian criticality with c = 1 which 
is separated from the sixfold-degenerate ordered and the 
disordered phases by two BKT-transition lines. Since 
these are found not to merge into a single curve, the in- 
termediate region as well as the self-dual line continues 
up to the critical decoupling point, and thus it can be 
regarded as a realization of the crossover phenomenon 
from the criticality c = I at the point (we shall discuss 
this issue in Sec. IIV|) |38j . In the previous researches, 
one expected that the critical region was terminated at 
a certain point, and a first-order phase transition be- 
tween the ordered and disordered phases occurred near 
the critical decoupling point [2ll [23, |32| • But, this possi- 
bility is removed. Instead, it is clarified that whereas the 
obtained phase diagram is, of course, quantitatively dif- 
ferent from that in the isotropic case, its structure, the 
stabilized phases, and the mechanisms of phase transi- 
tions are identical to those of the isotropic model 0, H3 • 
Therefore, we conclude that the present model belongs to 
the 6SC universality class, independently of the spatial 
anisotropy, and thus KO's assertion is confirmed. 

IV. DISCUSSIONS 

First, we discuss the nature of our phase diagram (in 
the isotropic case Qian and Blote performed the sim- 
ilar analysis The critical decoupling point with 
c = | becomes unstable against relevant competing per- 
turbations and exhibits crossovers to the behaviors con- 
trolled by the critical fixed points with lower symmetries 
|S8j . While one of those perturbations is the energy- 
density operator of the Ising model with the dimension 
1, another one may be a product of the magnetization 
operators on two of three sublattices which has the di- 
mension 2 x i |40|. Therefore, the crossover exponent 
4> = (2 — j)/{2 — 1) = | = 1.75 can predict the shape of 
the boundary around the point. In Fig.0thc log-log plot 
of the phase boundary line vi(u) is given by the triangles 
with the least-squares-fitting line to the data for the four 
smallest \K\\ [we define K% = ln(l + \/2)]. The inverse of 
the slope estimates the exponent </>f — 1.72. We also an- 
alyze v & &(u), i.e., the crossover line to the Gaussian fixed 
point with K = g (see the diamonds with the fitting line 




-5 -2.5 

ln(|K,|) 

FIG. 5: Log- log plots of the phase boundary line fi(w) and 
the self-dual line v a d(u) around the critical decoupling point 
[we define K% = ln(l + v2)]- The least-squares-fitting solid 
line for triangles (diamonds) estimates the crossover exponent 
<t> ¥ ^ 1.72 (4>af =± 1.76) in the F (AF) case. 

in the same figure). Then, the estimated exponent also 
takes the close value, i.e., 4>af — 1.76. Therefore, we 
can confirm the above crossover argument. On the other 
hand, the analysis of the BKT-transition lines Utj,l(w) 
becomes problematic. This may be mainly due to the 
existence of the corrections stemming from the marginal 
operators, while those are absent in the F case. 

Second, based on our results obtained in the above, 
we shall provide some comments on previous work: 
Miyashita, Kitatani, and Kanada performed the MC 
simulations at K 2 /Ki = 0.2 and 0.5 Hj. While they 
used some methods to estimate the transition points, 
the most specific ones are those by the MC-RG method 
|4ll |. They employed the real-space renormalization for 
each elementary triangle with corner spins Si, S2, S3, i.e., 
si + S2 + S3 — S1S2S3 1— ► s' 01 1 an d then evaluated 
two BKT-transition points. The comparison between 
Fig. |21 and their MC-RG data exhibits that the lower- 
temperature transition point at K2/K1 = 0.5 deviates 
from our phase boundary. One may attribute the dis- 
crepancy to the statistical errors in MC simulations, but 
there might be a possibility that the anisotropy effect 
could increase an uncontrollability of the real-space RG 
treatment. 

Pajersky and Surda employed the cluster TM method 
(a combination of the TM and a mean-field approxima- 
tion), and claimed that the intermediate critical phase 
possessed an incommensurate structure [22^ . They also 
implied that the lower-temperature transition between an 
incommensurate liquid and the commensurate ordered 
phase was the Pokrovsky-Tarapov type |43|. However, 
as exhibited, the dual sine-Gordon field theory @ well 
describes the phase transitions observed in this model, 
and there exists no relevant term to stabilize the incom- 
mensurate phase |44j . Further, since the duality relation 
is possessed by the effective model (u 7^ 0) and it inter- 
changes the transition points in upper and lower temper- 
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atures, they should be the same type. Therefore, their 
observations may be due to an artifact in their mean-field 
approximation treatment. 

de Queiroz and Dommany investigated the same model 
by the TM calculations [23j. They provided limited ev- 
idence for the existence of the BKT phase, and also 
alternative scenarios including an absence of the BKT 
phase. Meanwhile, based on the PRG calculation data 
and the scaled gaps, they drew the qualitative phase di- 
agram including the multicritical point and the direct 
transition line between ordered and disordered phases. 
However, it has been pointed out by several authors that 
the PRG calculations fail to estimate the BKT-transition 
points ^3 ; their results suffer from an inadequacy of the 
method. Further, they tried to check the KT criterion on 
their phase boundaries. They indeed estimated the expo- 
nent r\ (= 2xs) by the use of Eq. I|10(l and the fixed value 
C = 2/V3. But, as we have seen, £ depends upon model 
parameters, so the scaled gaps cannot give their univer- 
sal amplitude in the anisotropic case. Furthermore, their 
observation that the ratio of scaled gaps gives the flat- 
tening region of the parameter is understandable from 
our viewpoint that approximate cancellation of the geo- 
metric factor occurs by taking their ratio. Consequently, 
these exhibit that the nonuniversal quantity £ is impor- 
tant in the quantitative description using the numerical 
calculations, although it is not theoretically important. 

To summarize, we investigated the Kitatani-Oguchi 
model by the level-spectroscopy method in order to clar- 
ify spatial anisotropy effects and the global phase dia- 



gram. By taking into account the parameter dependence 
of the geometric factor properly, we analyzed the scal- 
ing dimensions of operators around the BKT-transition 
lines. Then, we numerically determined the phase dia- 
gram, where two types of Berezinskii-Kosterlitz-Thouless 
transition lines separate the intermediate critical phase 
from the ordered and disordered phases. Further, we 
evaluated the central charge and performed consistency 
checks among scaling dimensions in order to provide evi- 
dence of the universality class, and then we confirmed the 
assertion made by Kitatani and Oguchi. Some comments 
and comparisons with previous work were also given on 
the basis of our viewpoint and numerical results. 

While our approach has its basis on the argument of 
the Tomonaga-Luttingcr liquid observed in the ID quan- 
tum systems, it is widely applicable to the 2D classical 
systems with spatial anisotropy. 
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